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Abstract 

We examine how accurately the general HZV couplings, with V = Z, 7, 
may be determined by studying e + e~ — ► Hff processes at future e + e~ lin- 
ear colliders. By using the optimal-observable method, which makes use of 
all available experimental information, we find out which combinations of the 
various HZV coupling terms may be constrained most efficiently with high 
luminosity. We also assess the benefits of measuring the tau-lepton helici- 
ties, identifying the bottom-hadron charges, polarizing the electron beam and 
running at two different collider energies. The HZZ couplings are generally 
found to be well constrained, even without these options, while the HZ'y cou- 
plings are not. The constraints on the latter may be significantly improved 
by beam polarization. 
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1 Introduction 



The standard model (SM) of elementary-particle physics predicts a neutral scalar 
Higgs boson H as a remnant of the spontaneous breaking of its gauge symmetry. 
This particle is the only undiscovered ingredient of the SM so far. The experiments 
at the CERN Large Electron-Positron Collider (LEP2) were able to place lower 
bounds on its mass in the range 91.0-98.8 GeV at the 95% confidence level (CL) [|TJ. 
The search for the Higgs boson is a prime target of future colliders. Once the Higgs 
boson is found, its properties and interactions with other particles may be studied 
in detail with e + e~ linear colliders. If the Higgs boson is light, the bremsstrahlung 
process e + e~ — > HZ is expected to be the most promising process to study its 
properties and interactions and to search for deviations from the SM predictions. 

The purpose of this paper is to study systematically the sensitivities to general, 
non-standard couplings among the Higgs boson, the Z boson and a neutral vector 
boson V (V = Z,j). Since the Z boson has spin one, we take into account the 
angular distributions of its subsequent decays to fermion-antifermion pairs, in order 
not to loose information on the interference between amplitudes with different Z- 
boson helicities. On the other hand, we treat the Higgs boson as a final-state particle 
because it has spin zero. Thus, we study the production and decay processes 

e+e- - HZ- Z^ ff 

to obtain sensitivity to general HZV couplings. 

We first review previous studies on related problems. The angular distribution 
of e + e~ — > Hff has been analyzed for the SM at the tree level in 0. Expressions 
for the cross sections have been elaborated for beam polarization in [||. Radiative 
corrections have been investigated in [|J . A comprehensive review of the Higgs boson 
properties has been given in || . The HZZ form factors have been introduced in the 
study of composite light Higgs bosons M . Effects of the non-standard couplings have 
been discussed in 0, ||. Z-boson decay angular distribution in the process e + e" — > 
HZ have been analyzed as a means of distinguishing a scalar form a pseudoscalar 
Higgs boson in 0. 

We employ the optimal-observable method [[II], [XT], [T2], |T3[ to obtain constraints 



on the HZV couplings. This method provides the most efficient way to extract 
physical parameters from experimental data in the sense that the statistical errors 
on these parameters are minimized. Atwood and Soni have introduced optimal 
observable quantities in their analysis of electromagnetic form factors of the top 
quark |T0fl . The optimal-observable method has also been used in the measurement 



of the tau polarization [11]. This method has then been extended to the many- 



parameter case. It has been applied to the determination of the electroweak triple- 



gauge-boson fll2[ , Hti and HZZ couplings |L3 . 

This paper is organized as follows. In the next section, we cast the differen- 
tial cross sections of the process e + e~ — > Hff into a compact form, to which the 
optimal-observable method can be applied. We then discuss the properties of the 
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various terms therein under discrete symmetries. In Sect. |3|, we introduce an ef- 
fective Lagrangian for the HZV interactions and calculate the helicity amplitudes 
of the process e + e~ — > HZ. In Sect. |], we determine the achievable errors on the 
optimal observables introduced in the general expansion of the cross section. In 
Sect. ||, the optimal constraints on the effective coupling constants are discussed for 
typical experimental situations. Our conclusions are presented in Sect. |[ 



2 Cross section of e + e — > Hff 

In this section, we present the general angular distributions of the differential cross 
section of the production and decay process, 

P e i |) + e + (p s , -|) - Z*h\q) - H(p H ) + Z(p z , A), (2.1-a) 
Z(Pz, ^ - / (p„ y) + / (pj, -y) , (2-1-b) 

in a compact form suitable for the optimal-observable method. The four-momentum 
and helicity of each particle is indicated in parentheses; we have a = ±1, a' = ±1 
and A = 0, ±1 . 

We evaluate the production process ( |2.1-a|) in the centre-of-mass (CM) frame of 
the colliding beams. The production amplitude is then a function of the scatter- 
ing angle enclosed between the incoming electron and outgoing Z-boson three- 
momenta, p e and pz, respectively. The Z-boson helicity A is defined in the CM 
frame of the colliding beams. The y axis is chosen along the p e x pz direction. 
The decay process ( |2.1-b| ) is described in the rest frame of the outgoing Z boson. 
Here, the z-axis is chosen along the direction of pz (before the boost). The decay 
amplitude is a function of the polar angle 9 and the azimuthal angle ip of the / 
three-momentum. 

2.1 Angular distributions 

The angular distributions of the differential cross section of e + e~ — ► Hff may be 
written as 

. "J r , = t, [<r^(e,»,*>) + <^(e,M] • (2-2) 

a cos k3cl cos Udtp ~[ 1 J 

where c^ ,A ^ are mo del- dependent coefficients and F^ V ' A ^ are known functions of the 
angles G, 6 and (p. We shall present the definitions of c[ V ' A ^ and F^'^ below. The 
functions F^ V ' A ^ depend on the flavor of the final-state fermion / and the polariza- 
tion P of the initial-state electron; we have P = ±1 if the electron beam is purely 
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right/left-handed. We assume that the positron beam is unpolarized. In the deriva- 
tion of c\ V,A> and F^ V ' A \ we use the narrow- width approximation for the Z-boson 
propagator and the SM amplitude for the decay process Z — > //. 

We define reduced helicity amplitudes by extracting the angular dependence 
from the helicity amplitudes for e + e~ — > HZ as 



M*(e + e 



HZ) = MX A (6), 



(2.3) 



where 



<a=o( ) = ~^ as[ne > 



<A= ± (e) 



;(1 + o-Acos0). 



(2.4) 



The amplitudes M^ do not depend on 0. 

The coefficients c\ V ^ and cf^ are expressed in terms of the amplitudes M^ as 
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Here, the + (— ) sign refers to V (A), and the subscript R (L) stands for a = +1 
(a = —1). These eighteen coefficients contain all observable consequences of the 
reduced amplitudes. 

The functions are defined as 



(v 



= - sin 2 sin 2 9, 
4 

= — (l + cos 2 0)(l + cos 2 #) -^-COS0COS#, 

16 4 

r rPAf 

= sin 20 sin 26 cos p H — L sin sin 6 cos ip, 

16 4 

= - sin 2 sin 2 6 cos 2(p, 

r rP 'At 
= ——sin 20 sin 26 sin <p H — sin sin 6 sin ip, 



(2.6-a) 

(2.6-b) 

(2.6-c) 
(2.6-d) 

(2.6-e) 
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P 6 (y) = ^ sin 2 sin 2 # sin 2^, 
8 



F (V) = 



Pi^ = sin 26 sin 9 cos ip — — sin 6 sin 29 cos ip, 



-^(1 + cos 2 0) cos 9 + — cos 0(1 + cos 2 9), 
8 8 

rA f . . „ rP 



f rP 
- sin 20 sin 9 sin </? sin sin 26 1 sin ip. 



(2.6-f) 
(2.6-g) 
(2.6-h) 
(2.6-i) 



The common coefficient r contains some phase space factors and the branching 
fraction of the Z — > // decay, 



r = - 



1 fez 3 
4 32vrs 4tt 



and Ay is the left-right asymmetry of this decay, 



2 • 



(2.7) 



(2.8) 



Here, s is the square of the CM energy, and (3hz is the two-body phase-space factor, 



P. 



HZ 



\ 



(2.9) 



The functions p/^ are obtained by flipping the P dependence. Specifically, if we 



write P 



(V) 



ai + /3«P, with P-independent functions ai and fa, then we have 



ctjP + fa. The explicit expressions read 
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The three angular functions Fi A \ and F^' vanish if P = 0. One could measure 



(A) 



and c^' if \P\ ^ by combining experiments with opposite polarizations 



1A) 



P = \P\ and P 



I PL 



For most of the hadronic decay modes of the Z boson, the final-state fermions 
/ and / cannot be distinguished. Then, we have to average over the configurations 
with (B, 9, ip) and (0, 7r — 9, <p ± it) as 



<V,A), 



(e,9,<p) + F l {v > A \&,7r-9,<p±ir) 



(2-11) 



This corresponds to setting Af = if / is a quark. If P = 0, then one can measure 
the coefficients cf^, . . . , and Cj A \ . . . , cj, , while the coefficients c^, . . . , Cg 
and c^, . . . , c 6 are only measurable if P ^ 0. 

When the Z bosons decay to neutrino pairs, one can only measure the distri- 
bution, so that the (9, tp) dependences should be integrated out. Then, we have 

da 
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(2.12) 

(2.13-a) 
(2.13-b) 

(2.13-c) 

(2. 13-d) 

(2.13-e) 

(2.13-f) 



while the other F 



(V,A) 



functions vanish. 



2.2 Discrete symmetries 

We now discuss the properties of the F^ V,A ^ functions under the discrete symmetries 
CP and CPT. Here, T is the naive time reversal symmetry, which flips the momen- 
tum and spin of all the particles, but does not reverse the time flow from the initial 
state to the final state. The non-vanishing of the CPT-odd coefficients is related to 



the presence of absorptive parts in the amplitudes fL4 |. 

The electron beam polarization P, the decay asymmetry Af and the angular 
variables transform under the discrete symmetries as 

(P, Af, 0, 9, if) ™ (P, Af, vr - 0, 7T - 9, 2ir-<p), (2.14-a) 

(P, A f ; 0, 9, ip) t (P, Af 0, 9, 2tt-<p). (2.14-b) 
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We can then obtain the symmetry properties of the ' functions, which are 
summarized in Table |I|. 



Table 1: CP and CPT properties of the F± ' functions. A + (— ) sign means even 
(odd) under the symmetry. 



i 


1 


2 


3 


4 5 6 


7 8 9 


CP 


+ 


+ 


+ 


+ - - 


- - + 


CPT 


+ 


+ 


+ 


+ + + 





The coefficients have the same symmetry properties as the functions F^ V,A \ 
The coefficients and c^'^ are sensitive to CP-odd and CPT-even quantities. 

When the new-physics effects are generated by exchanges of heavy particles, then 
the induced vertices should be CPT even. The three CPT-odd coefficients , 
Cg and Cg^'^ should be proportional to the absorptive parts of the amplitudes 
which contain light particles in the loops. 

3 Helicity amplitudes for e + e~ — ► HZ 

In this section, we first introduce general couplings and effective form factors for the 
HZZ and HZ'j interactions. We then present the helicity amplitudes of e + e~ — > HZ 
using these form factors. 

We adopt the effective HZV interaction Lagrangian from ||. It reads 



C ea = (l + a z ) 9 -^HZ^ + ^- £ \b v HZ^ 

+ c v {d^HZ v - dvHZ^) + byHZ^V^] , (3.1) 



where V^ u = d fl V u — d u V^ and V^ u = e^apV" 13 with the convention e i23 = +1- We 
have neglected the scalar component of the vector bosons by putting 

= d^V = 0. (3.2) 

Then, the most general parameterization of the HZV interaction involves seven 
couplings, az, bz, cz, b 7 , c 7 , bz and 6 7 , which are constants as long as we only 
consider operators through mass dimension five. We note in particular that the 
operator identity 

HZ^Z^ = HZ^Z^ 



- l -HZ^ - l -{d,HZ v - d u HZ^)Z^ (3.3) 
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holds under the condition (|3T2|). The five couplings az, bz, cz, & 7 and c 7 are CP 
even, while the remaining two couplings, bz and 6 7 , are CP odd. In the effective 
Lagrangian ( |3.1| ), we have factored out the Z-boson coupling g z and appropriate 
powers of m z to render the couplings dimensionless. In the SM, we have az = by = 
cv = by = at the tree level. 

The form factors for the generic HZ a Vp vertex may then be written as 



T a0(<liPz) =9z m Z 



K [s)g a p + —q a p zi 3 + —tafS^Pz 



m z ' m z 



(3.4) 



where the virtual y-boson momentum q is taken to be incoming and the Z-boson 
momentum p z to be outgoing, as depicted in Figure [I] and process fl2.1|), and s = q 2 . 
All form factors are dimensionless functions of s. The four form factors /if , /if, 
hj and hi are CP even, while the two form factors /if and /ig are CP odd. It is 
straightforward to express the form factors in terms of the seven couplings of the 
effective Lagrangian (|3.1|) : 



/if (s) = 1 + a z ) + 2c z ^ + 2{b z - c z ) ^ ^, 3.5-a 

m z m z 

h z 2 {s) = -A(b z -c z ) } (3.5-b) 
/if (s) = -46 z , (3.5-c) 

/i7(s) = 2c 7 ^ + (6 7 - c 7 ) g + m| ~ m ^ ^ ( 3 - 5 " d ) 

/i](s) = -2(6 7 -c 7 ), (3.5-e) 
/i^(s) = -2& 7 . (3.5-f) 

Although the effective Lagrangian has seven couplings, there are only six form fac- 
tors. Thus, one combination of couplings cannot be measured at one given collider 
energy. Details will be discussed in Sect. |E| 

We now evaluate the helicity amplitudes M^(e + e~ — > HZ) for the production 
process. After extracting the angular dependence according to ([2. 3D, we obtain the 
reduced amplitudes^ as functions of s: 



M*=°(s) = -9z9:V2!EzDz(s) ( hf + hi^& 



+ eg z V2~sE z D^s) ( hj + hl ^f* ] , (3.6-a 



z 



M*=%s) = -g z g:V2Sm z D z (s) (/if + zA/if 



+ eg z V2lm z D 7 (s) (hj + i\hJ ^ E f Z ) . (3.6-b) 



1 In there is a misprint in the relative sign of the helicity amplitudes. There should be an 
overall minus sign on the right-hand side of (2.3) therein. Furthermore, Im [(M+ + M~)(M°)*] 
and Im [(M+ — M~)(M°)*] should be interchanged in the first column of Table 1. 
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where g e R = sin 2 9w, <?£ = — 1/2 + sin 2 6w, Pz = y 1 — m2 z/E z an d 

= 2 I ■ r , (3-7) 

D y (s) = -. (3.8) 
The energy of the outgoing Z boson is 

Et = £( 1 _ . (3 . 9) 



The couplings g>fj and g£ have almost the same magnitudes, but their signs are op- 
posite to each other. Thus, the coefficients are sensitive to the HZZ couplings, 
while the coefficients c A are sensitive to the HZ'y couplings. 



4 Optimal observables 



We now employ the optimal-observable method [12, |H| to obtain the errors on the 
coefficients q and c~ A \ In order to simplify the notation, let us re-sequence the 
coefficients and functions for the time being as 

( Cl ,...,c 18 ) = (c?\...,4 v \c[ A \..., C W), (4.1-a) 
(F X , . . • , F 18 ) = (Fi V \ . . . , *f\ f} A \ . . . , F^). (4.1-b) 

According to the optimal-observable method, the covariance matrix Vy for the co- 
efficients Cj is given by 

=L — d cos6 cos 6d<p, 4.2 

where L is the integrated luminosity of the experiment and 

Ssm(6, 6, tp) = d ° SM (9, 6, <p). (4.3) 
d cos 6 cos 



The statistical error on q is v^4i; an d the correlation between the error on q and 
that on Cj is V^-/ JVaVjj- We use this method to obtain V^ 1 for each decay mode 
of the Z boson. We then combine these results for all the Z-boson decay modes. 
The differential cross section £sm is a linear combination of F[ V \ . . . , ff > and 
, which are CP and CPT even. The component of the covariance 
matrix vanishes if Fi and Fj have different CP or CPT properties. Thus, the 
covariance matrix is block diagonalized into four sub-matrices according to the CP 



and CPT properties of the Fi functions discussed in Sect. 2.2. Notice that this 
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argument is only valid if we integrate in ( f4.2| ) over the full angle domains. In 
practice, there are excluded regions due to the geometry of the detectors or cuts for 
selecting events. Thus, the block diagonal structure of the covariance matrix is only 
approximately realized in practice. In the present study, we shall integrate over the 
full phase space. 

We estimate how the optimal errors are reduced by the following three additional 
techniques. The first one is the tau helicity measurement. We adopt e r = 40% as 
the efficiency factor to determine the helicities of the decaying r + or r~ leptons. 
The second technique is the electric-charge identification for the bottom quarks and 
antiquarks. The charge of a hadron B containing one b or b quark can be identified 
via the decay mode B — > lu+X . We assume an efficiency of = 20% for identifying 
the charges of the decaying b or b hadrons. The third technique is to employ electron 
beam polarization P. We take \P\ = 90% as the target polarization. Specifically, 
we assume that one half of the beam is polarized with P = 0.9 and the other half 
with P = -0.9. 

For the fraction e r of the Z — > t + t~ decays, one can distinguish the tau polariza- 
tion. In order to assess the possible benefits of the tau polarization measurement, we 
make the simple assumption that the efficiency for observing a left- or right-handed 
tau lepton is e T . Then, we substitute in ( |2.6|) and ( |2.10| ) 



3/3 



HZ 



512tt 2 s 
A f = ±l 



Br(Z 



T T 



9l/r 



{sir + (9 t r) 



(4.4-a) 
(4.4-b) 



for left/right-handed tau leptons. In actual experiments, one should not only es- 
timate the efficiency factor e T for each pair of r + and t~ decay modes, but also a 
correlation between the constraints from production and those from produc- 
tion. We return to this problem at the end of this section. 

Throughout our numerical analysis, we set m z = 91.187 GeV []15|| , a = 1/128.9 

sin 2 9 W = 0.2312 and g z = yA\f2G F m 2 z = 0.74070. As an example, we show 



results for the Higgs boson mass m# = 120 GeV, the CM energy y 7 ? = 250 GeV 
and the nominal integrated luminosity L = 10 fb _1 . 

We first present the results for e T = e& = P = 0. The results for the CP-even 
and CPT-even coefficients and their block in the covariance matrix are 
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The results for the CP-odd and C PT-even coefficients are 
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The results for the CP-odd and CPT-odd coefficients are 
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The results for the CP-even and CPT-odd coefficients are 
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c K 8 v ' and Cg^ are relatively large because 



and Pi" 4 '* vanish if 



There are no constraints on c\ ' . r 
P = 0. The errors on c| , C3 ', c 

the corresponding F^ V ' A ^ functions in ( |2.10| ) are suppressed by the smallness of 
and the vanishing of P. 

Next, we present the results for e T = 40%, e& = 20% and |P| = 90%. The results 
for the CP-even and C PT-even coefficients are 
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The results for the CP-odd and C PT-even coefficients are 



c iV) 
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± .0026 
0± .0030' 
0± .0029 
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The results for the CP-odd and CPT-odd coefficients are 

Oi.0011 / 1 
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.17 
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The results for the CP-even and CPT-odd coefficients are 

<f°=0±.0023 / 1 \ 

4 A) =0±.002l' {-.13 1 J' (412) 

The errors on c A with i = 1, . . . , 6 and q with i = 7, 8, 9 are reduced to the 
level of those on the other coefficients because the suppression of the corresponding 
p(Y> A ) f unc tions in ( |2.10| ) is weaker than in the case of e T = e& = P = 0. On the 
other hand, the errors on the other coefficients are not further reduced relative to 
the former situation. 

We mention here the effect of the correlation between the constraints from t r 
production and those from r L production. In actual experiments, t l and tr leptons 
can only be identified on a statistical basis. The analyzing power of the semileptonic 
tau decays is, in principle, equal for all the semileptonic tau decay modes [0. By 
using the r~ — > v T Tt~ decay mode, we evaluate the effect of the t l -t r correlation 
on the errors on c\ V ' A \ We find that the errors on c\ V,A ^ may be increased by about 
20% in actual experiments. 



5 Constraints on general HZV couplings 

We are now ready to study the sensitivities to the seven general HZV coupling 
constants. The errors on these couplings are obtained from those on c\ V,A ^ by using 
( f2~5"D , (|3.5|) and ( |3.6|) . We quantitatively analyze the usefulness of electron beam 
polarization and of an additional experiment with another beam energy. For con- 
sistency of the analysis that includes the operators through mass dimension five, we 
only keep in c[ V ' A ^ terms linear in the coupling. 



5.1 Real part 

We first discuss the constraints on the real parts of the general HZV couplings. 
The constraints on the real parts of the CP-even couplings are obtained from the 
CP-even and CPT-even coefficients c{' A \ . . . , c± , while those on the real parts 
of the CP-odd couplings are obtained from the CP-odd and CPT-even coefficients 

c^ ,A) and c§' A \ 

We first present the results for ^fs = 250 GeV. The optimal errors on the HZV 
couplings are summarized in Table [| We only gain sensitivity to six combinations 
of couplings. As long as we consider experiments at a fixed collider energy, one 
combination of couplings cannot be measured. The unmeasurable combination of 
couplings is determined from ( p.5|) and reads 

2 

TYl 
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Table 2: Optimal errors on the real parts of the general HZV couplings at y/s = 
250 GeV. 







0.4 






0.4 


e& 






0.2 




0.2 


|P| 








0.9 


0.9 


Re{b z + .059a z ) 


.0061 


.0036 


.0033 


.0030 


.0029 


Re (c z + .059a z ) 


.013 


.0076 


.0070 


.0061 


.0061 


Re 6 7 


.19 


.072 


.053 


.0085 


.0084 


Re c 7 


.12 


.047 


.035 


.0053 


.0052 


Re ~b z 


.012 


.011 


.010 


.010 


.0091 


Re 5 7 


.094 


.036 


.026 


.016 


.013 



It is independent of the final-state fermion flavour / and the electron beam polariza- 
tion P. We are thus insensitive to this combination for all Z-bosons decay modes. 
Since a z is the dominant part of ( |5.1| ), we fix a z to obtain the optimal sensitivities 
to the remaining six coupling constants b z , c z , 6 7 , c 7 , b z and b 1 . The combinations 
Re (bz + .059az) and Re (c z + .059a^), which appear in Table |2], are orthogonal to 
the unmeasurable combination fl5.1|) . 

For e T = £{, = P = 0, we have good sensitivities only to the three HZZ couplings 
bz, cz and bz, but not to the HZj couplings 6 7 , c 7 and 6 7 , which is evident from 
Table || The functions F[ ' with % = 1, . . . , 6 are suppressed in magnitude by the 
smallness of Af and the vanishing of P, while the functions F[ v ^ with i — 1, ... ,6 
have unsuppressed parts. 

By using any of the three additional techniques, we gain better sensitivities to 
the HZ'y couplings because F^ with % = 1, . . . , 6 are then less suppressed. The 
measurement of the tau helicity with 40% efficiency reduces the errors on 6 7 , c 7 and 
b-y by a factor of about 2/5 relative to the case without tau helicity measurement. 
Bottom charge identification with 20% efficiency reduces the errors on these cou- 
plings by a factor of 2/7. We observe that the tau helicity measurement leads to an 
improvement comparable to that for the bottom charge identification. This may be 
understood qualitatively from the relation 



e b Br(Z -►&&).„. 2/5 , . 

^ 53F^) '*' ~ L5 ~ W (5 - 2) 

The electron beam polarization is the most efficient technique for improving the 
sensitivities. It reduces the errors on the CP-even (CP-odd) HZ'y couplings by a 
factor of about 1/20 (1/6). A qualitative understanding hereof is obtained from the 
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relation 



1 



|P| « 8.2 



2/5 



(5.3) 



Br(Z -> r-r+) r 1 " " 1/20' 

for the CP-even couplings. We find from Table ^| that the errors on the CP-even 
HZZ couplings are reduced by a factor of 1/2 with these three additional techniques, 
while the CP-odd HZZ couplings are almost unchanged. 

For e T = = P = 0, the errors on the real parts of the couplings and the 
corresponding correlation matrix are 



Re (b z + .059a z ) 
Re (c z + .059a z ) 

Re 6 7 

Rec 7 

Reb z 

Re b-y 



0± .0061 
0± .013 
0± .19 
± .12 ' 
0± .012 
± .094 



/ 1 



\ 



.95 1 



1 



.83 -.89 -.99 1 

V 



(5.4) 



-.46 1 / 



There are strong correlations among the errors on bz, c z , b 1 and c 7 . Thus, a cer- 
tain combination of parameters is more strongly constrained than the individual 
parameters. The eigenvector with the smallest eigenvalue and its error read 



Re (.08a z + .906 z + .42c z + .056 7 + .08c 7 ) = ± .00076. 



(5.5) 



The above combination may be understood qualitatively by observing that an exper- 
iment at y/s = 250 GeV operates near the threshold of the HZ production process, 
where (3hz ~ 0. Near the threshold, the form factors /if and K[ play a dominant 
role in the helicity amplitudes, while the residual form factors are suppressed by the 
smallness of Phz- Specifically, we have 



Mi 



g z V2sm z -g z g*D z (s)hf + eD^hJ + 0(J3) 



(5.6) 



Furthermore, the sensitivities to the HZ'-f couplings are diminished for e T = e& = 
P = because sin 2 (9^ w 1/4, Thus, near threshold there is good sensitivity to the 
following combination of couplings: 



hf w 1 + a z + 4b z 1 + 



rrir 



+ 2c z - 



m 



Tflr 



(5.7) 



This is essentially the combination that appears in the constraint ( |5.5| ). 

In models with multiple Higgs doublets, including the minimal supersymmetric 
extension of the SM (MSSM), the coupling az is modified at the tree level, while 
the couplings by, cy and by only receive corrections at the loop level. Thus, we 
discuss here the sensitivity to az when by = cy = by = 0. From (|5.4|) , we obtain at 
v/i = 250 GeV 

a z = ±0.010 (5.8) 

if by = Cy = 0. 
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For e T = 
found to be 



40%, e& = 20% and \P\ — 90%, the errors and correlation matrix are 



Re (b z + .059a z ) 
Re(c z + .059a z ) 

Re 6 7 

Rec 7 

Re~b z 

Re 6 7 



± .0029 
± .0061 
± .0084 
0± .0052' 
± .0091 
0± .013 



/ 1 
-.96 1 
-.08 .08 1 
.08 -.09 -.99 1 

V 



(5.9) 



-.09 1 / 



We obtain similar correlation matrices for the other situations when only one of 
the three additional measurements is employed. The correlations between the HZZ 
and HZ'-y couplings then disappear. There are still strong correlations between b z 
and c z and between 5 7 and c 7 . The eigenvectors of the two smallest eigenvalues and 
their errors are 

Re (.025a z + .2Sb z + .14c z + .506 7 + .81c 7 ) = ± .00068, (5.10-a) 
Re (.074a z + .866 z + .40c z - .176 7 - .27c 7 ) = ± .00078. (5.10-b) 

As in (|5.6| ), the form factors hf and h\ become dominant near the threshold, and 
there is good sensitivity to the following two combinations of couplings: 

hf ^l + a z + 4b z (l + — ) + 2c z ^f, (5.11-a) 
\ m z ) m z 

h^2bJl + ^)+2cJ^ + r ^). (5.11-b) 



rrir 



rrir 



rrir 



The most-strongly constrained combinations listed in (|5.10|) are essentially linear 
superpositions of hf — 1 and hj as given in ( |5.11| ). 

We now discuss the sensitivity to a z when by = cy = by = 0. The six-parameter 
constraints ( |5TP| ) then lead to the one-parameter constraint 



a z = 0± 0.010 



(5.12) 



at a/s = 250 GeV. This constraint is same as in ( |5.8|) . The error on a z is not 
diminished by any of the three experimental options. 

Figure || displays the contours of \ 2 — 1 (39% CL) in the (6 7 , c 7 ) plane for 
the different modes of experiment. The other five degrees of couplings have been 
integrated out. We observe that there is a strong correlation between 6 7 and c 7 
for all experimental methods. As mentioned above, the specific combination of b 1 
and c 7 contained in ( |5.10| ) is thus tightly restricted. We can see from Figure ^ that 
the individual sensitivities to 6 7 and c 7 are drastically improved by the electron 
beam polarization. This means that we can obtain strict constrains on any model 
that predicts large HZ'-f couplings by using data from experiments with polarized 
electron beams. 



15 



Figure [3] displays the contours of % 2 = 1 in the {b z , c z ) plane. We see that the 
HZZ couplings are well constrained even if e T = % = P = 0. The three charge and 
polarization measurements lead to moderate reductions of the errors on b z and c z . 

Figure ^ shows the contours of \ 2 — 1 m the (6^,6 7 ) plane. The three charge 
and polarization measurements mainly reduce the error on b 1 . The reduction of the 
error on b 1 is transferred to that on bz via the correlation between bz and b 1 . 

Next, we consider the case of yfs = 500 GeV. The results for e T = 40%, e& = 20% 



and P = 90% are 




Re{b z + M6a z ) 


= 0±.0015 


Re(c z + MQa z ) 


= ± .0007 


Re 6 7 


= ± .0024 


Rec 7 


= 0±.0005' 


Reb z 


= ± .0042 


Re 6 7 


= ± .0052 


Although the cross section at \fs 



( 1 

-.77 1 
-.09 .07 1 
.07 -.09 -.84 1 
1 
V -.09 1 J 



(5.13) 



500 GeV is smaller the one at y/s = 250 GeV, 
the errors are reduced because b v , c v and by are accompanied by a factor of s/m 2 z . 
Increasing the CM energy from = 250 GeV to ^/s = 500 GeV reduces the errors 
on bz and b 1 by a factor of 1/2 to 1/3 and those on cz and c 7 by a factor of 1/7 to 
1/10. This also reduces the errors on the CP-odd couplings bz and 6 7 by a factor 
of 1/2. The strong correlations are lost because the experiment at y/s = 500 GeV is 
far above the threshold, and the various form factors in the helicity amplitudes are 
non- negligible. At \/s = 500 GeV, the sensitivity to az becomes 



a z = ±0.021 



(5.14) 



if by, cy and by are fixed to zero. This error is larger than the one at yfs = 250 GeV 
because a z is the coefficient of the renormalizable dimension-four operator. We 
conclude that az may be well measured at the CM energy where the cross section 
of HZ production has its maximum. 

Finally, we present optimal constraints on the seven parameters in ( |3.1|) by com- 
bining the analyses at ys = 250 GeV and = 500 GeV with e T = 40%, e b = 20% 
and P = 90%. Here, we face the problem that our results depend on the integrated 
luminosities at the two energies. At fixed energy, our constraints scale as L -1 / 2 . We 
can arbitrarily scale our results by changing the nominal value of L, which could 
include more realistic experimental efficiencies. Once we combine the analyses at 
the two energies, our results will depend on the ratio of the two respective values 
of L, which we cannot fix a priori. For simplicity, we assume the same luminosity, 
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L = 10 fb 1 , at both energies. The results are then 



Re{b z + .066a z ) 
Re cz 
Re6 7 
Rec 7 
Reb z 
Re 6^ 



± .0009 
± .0006 
0± .0015 
0± .0004' 
± .0038 
± .0049 



/ 1 
-.68 1 
-.08 .07 1 
.06 -.08 -.79 1 

V 



1 

-.09 1 J 



(5.15) 



The minimum x is found to be 



X 2 - 

Amin 



az V 
0.024/ 



(5.16) 



Re az 
Reb z 
Recz 
Re b. 
Re a. 



7 



■"T 

Re6 z 
Re b\ 



± .024 
0± .0018 
± .0006 
0± .0015, 
0± .0004' 
± .0038 
± .0049 



So far, we have considered az as a fixed parameter. Now, Xmm * s a function of az, 
so that we can obtain optimal constraints on all seven HZV couplings. The result 
is 

/ 1 
-.87 1 
-.02 -.31 1 
.00 -.04 .07 1 
.00 .03 -.08 -.79 





V 



1 






1 



.09 1 



(5.17) 



Now, az is weakly constrained. There is a strong correlation between az and bz- 
This reflects the fact that the combination Re (bz + -06Qaz) in ( 5.15| ) is more strongly 
constrained than Re bz ■ 

Figure [5] illustrates how a z is constrained. As mentioned above, there is a combi- 
nation of couplings, namely the one in ( |5.1|) , that is not constrained by an experiment 
at a single CM energy. The projections of the cylinder defined by x 2 = 1 onto the 
(az, bz), (bz, cz) and (cz, az) planes are indicated as the stripes between the dashed 
(thin solid) lines for yfs = 250 (500) GeV. Because the direction of the cylinder 
varies with ^/s, the measurements at the two energies, a/s = 250 GeV and 500 GeV, 
lead to individual constraints on az, bz and cz- 

So far, we have assumed that az is constant, 
energy dependence, 

a z (s) = a z (0) + sa' z (0) -\ 



In general, az may have some 



Oi 



(5.18) 

The 0(s 2 ) term is neglected in our approximation. In terms of operators, the deriva- 
tive term corresponds to a dimension-five operator and should be exhausted by the 
effective Lagrangian ( |3.1| ). In fact, it may be written as a linear combination of the 
b z and Cz terms because of the operator identity (|3T 
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5.2 Imaginary part 



We now discuss the sensitivities to the imaginary parts of the general HZV cou- 
plings. Imaginary parts can arise if the couplings are induced by new interactions 
involving particles that can be produced at the energy of the considered experiment. 
For consistency, we neglect the absorptive part in the Z-boson propagator as we ig- 
nore all absorptive parts in the SM amplitudes. The constraints on the imaginary 
parts of the CP-even couplings are then obtained from the CP-even and CPT-odd 
coefficient c^' A \ while the constraints on the imaginary parts of the CP-odd cou- 
plings are obtained from the CP-odd and CPT-odd coefficients C7 and c§ . 
The results are summarized in Table |3|. In contrast to the real parts of the HZV 



Table 3: Optimal errors on the imaginary parts of the general HZV couplings at 
y/s = 250 GeV. 







0.4 






0.4 








0.2 




0.2 


|P| 








0.9 


0.9 


Im (b z - c z ) 


.25 


.095 


.071 


.015 


.014 


Im (bj — c 7 ) 


.055 


.027 


.023 


.018 


.018 


Im bz 


.049 


.018 


.014 


.0026 


.0026 


Im6 7 


.010 


.0050 


.0043 


.0032 


.0032 



couplings, one can only measure four combinations of their imaginary parts. The 
other combinations, Im az, Im (bz + Cz) and Im (fe 7 + c 7 ), only affect the form factor 
h\ , but not 1%2 or h\[ . They contribute to the amplitudes as a common overall 
phase, and hence they do not alter the extracted values of c\ V ' A \ Thus, these other 
combinations cannot be measured. 

For e T = €b = P = 0, the functions f[ V ' with i = 7, 8, 9 are suppressed in 
magnitude by the smallness of Af and the vanishing of P, while the functions 
with i — 7, 8, 9 have unsuppressed parts. Thus, the errors on Im (6 7 — c 7 ) and Im6 7 
are much smaller than those on Im (bz ~ cz) and Imfe^, respectively, as may be seen 
in Table § 

By using any of the three charge and polarization measurements, we gain better 
sensitivities to the HZZ couplings because the functions P^ with i = 7, 8, 9 are 
less strongly suppressed. The measurement of the tau helicity with 40% efficiency 
reduces the errors on the HZZ couplings by a factor of about 2/5. The bottom 
charge identification with 20% efficiency leads to a reduction by a factor of 2/7. 
The electron beam polarization is the most efficient technique for improving the 
sensitivities. It reduces these errors by a factor of 1/20. The errors on the HZ'y 
couplings are reduced by a factor of 1/2 with tau helicity measurements or bottom 
charge identification, and by a factor of 1/3 with beam polarization. 
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For e T = e& = P = 0, the errors and the correlation matrix are 



Im {b z - 
Im (bry — 
Im b z 
Im b\ 



0± .25 
0± .055 
0± .049' 
0± .010 



( 1 
-.94 1 



V 



(5.19) 



-.95 1 / 



There are strong correlations between the errors on the first two terms and those 
of the latter two. The eigenvectors of the two smallest eigenvalues and their errors 
read 



.20 Imbz + .98 Im& 7 = ± .0031, 

.20 Im (b z - c z ) + .98 Im (6 7 - c 7 ) = ± .018. 



(5.20-a) 
(5.20-b) 



For e T = 40%, e b = 20% and P = 90%, we have 



Im (b z -c z ) = 0± .014 
Im(6 7 -c 7 ) =0±.018 
lmb z =0±.0026' 



Im b~ 



0±.0032 



( 1 
-.10 1 



V o o 



-.10 1 / 



The strong correlation between the HZZ and HZ'-f couplings is lost for 
charge and polarization measurements. 

Next, we consider the CM energy y/s = 500 GeV. For e T = 40%, e& = 
\p\ = 90%, we find 



(5.21) 

the three 
20% and 



Im (b z - c z ) = ± .0033 
Im (6 7 - c 7 ) = ± .0037 
Im6 z =0±.0015' 



Im b*, 



0±.0017 



( 1 
-.10 1 



V o o 



\ 



1 

-.10 1 / 



(5.22) 



Although the cross section at y/s = 500 GeV is smaller than the one at y/s = 
250 GeV, the errors on the couplings are reduced because of the s/m% factors 
multiplying the above couplings. The three combinations Ima z , Im (bz + cz) and 
Im (b y + c 7 ) cannot be measured even if the CM energy is varied. 



6 Conclusion 

In the present paper, we have performed a systematic study of the angular distri- 
butions of the process e + e~ — > Hff in order to assess the sensitivities to the seven 



general HZV couplings by using the optimal-observable method [H], [12], 0. To 
that end, we have expanded the differential cross section as a sum of the products 
of the eighteen model-dependent coefficients c^'^ , which contain all the dynamical 
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information on the HZV couplings, and the corresponding eighteen angular func- 
tions F^ V,A \ which depend on the production and decay kinematics, the final-state 
fermion flavor /, the tau polarization and the electron beam polarization P. 

As for the real parts of the HZV couplings, one can only measure six combi- 
nations at a given CM energy y/s. At y/s = 250 GeV, we gain optimal errors of 
order 1 x 10~ 2 (1 x 10 _1 ) for the HZZ (HZ~/) couplings assuming L = 10 fb _1 
and m H = 120 GeV. A tau helicity measurement with 40% efficiency reduces the 
optimal errors on the HZ'-f couplings to about 2/5 of those obtainable without such 
a measurement. A bottom charge identification with 20% efficiency reduces these 
errors to about 2/7 of those with unidentified bottom charge. An electron beam 
polarization of 90% reduces the optimal errors on the CP-even (CP-odd) HZ'-f cou- 
plings to about 1/20 (1/6) of those with unpolarized beams. The reduction of the 
errors on the HZZ couplings is at most by 1/2. The sensitivities to the HZV cou- 
plings depend on y/s. The errors on the real parts of the HZV couplings decrease 
by a factor of about 1/2 to 1/10 when one increases y/s from 250 GeV to 500 GeV. 

As for the imaginary parts of the HZV couplings, we can only measure four 
combinations, as long as we only keep terms linear in the couplings. Without the 
three charge and polarization measurements, we achieve optimal errors of order 
1 x 10~ 2 for Im& 7 , 5 x 10 -2 for Im (6 7 — c 7 ) and lmb Z} and 3 x 10 _1 for Im (b z — cz) 
with L = 10 fb~ x at y/s = 250 GeV. The optimal errors on the HZZ couplings are 
reduced by factors of about 2/5, 2/7 and 1/20 with the tau helicity measurement, 
the bottom charge identification and the electron beam polarization, respectively. 
The errors on the HZ'-f couplings are at most diminished by a factor of 1/3 with the 
three charge and polarization measurements. When one increases y/s from 250 GeV 
to 500 GeV, the errors decrease by a factor of 1/2 to 1/5. 

In our analysis, we have considered the general HZZ and HZ'-f interactions. We 
have neglected the contribution from the dimension-five HZee operator, 

— £ g? Zee HZ»e llM P a e, (6.1) 

m Z <xt± 

with P± — (1 ± 75)/2, which contributes to the cross section at the same order 



as the operators in the effective Lagrangian ( |3.1| ). The simple and very general 
treatment of the observables in Sect. ^ is no longer valid if the terms in (|6.1| ) are 
significant. The optimal constraints on the effective couplings can still be obtained 
by directly studying the g^ Zee dependences of the differential cross sections. We 
believe, however, that our approach will be useful in constraining theories that affect 
the HZV couplings more significantly than the HZee couplings. The contributions 



from the third-generation squarks in the MSSM |L8[ provide one an example. 
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Figure captions 



Figure [l|: General HZV coupling. The arrows indicate the direction of the four- 
momentum flow. 

Figure [2|: Contours of x 2 = 1 in the (6 7 , c 7 ) plane for yfs = 250 GeV. The other 
degrees of general couplings are integrated out. The central values of 6 7 and 
c 7 are assumed to coincide with their SM values, 6 7 = c 7 = 0. The errors 
are estimated by means of the optimal-observable method under the following 
four conditions: (a) e T = e& = P = 0; (b) e r = 0.4 and e& = P = 0; (c) e& = 0.2 
and e T = P = 0; or (d) e r = 0.4, e b = 0.2 and \P\ = 0.9. 

Figure [3|: Contours of x 2 — 1 m the (bz, cz) plane for yfs = 250 GeV. = and 
the other degrees of general couplings are integrated out. The central values 
of bz and cz are assumed to coincide with their SM values, bz = cz = 0. 

Figure [|]: Contours of x 2 = 1 in the (bz, b 7 ) plane for ^fs = 250 GeV. The other 
degrees of general couplings are integrated out. The central values of bz and 
6 7 are assumed to coincide with their SM values, bz = b 7 = 0. 

Figure [5|: The projections of the \ 2 — 1 contours onto the (a z ,b z ), (bz-.Cz) and 
{cz, a z) planes for ^Js = 250 GeV and 500 GeV. The central values of a z , 
b z and cz are assumed to coincide with their SM values, a z = b z = c z = 0. 
The errors are estimated under the condition e r = 0.4, e& = 0.2 and \P\ = 0.9. 
Dashed (thin solid) lines represent the one-cx contours obtained at y/s = 250 
(500) GeV. Thick solid curves represent the combined one-cx contours. 
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Figure 2: 
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